Modern Michelson-Morley experiment using cryogenic optical resonators 
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We report on a new test of Lorentz invariance performed by comparing the resonance frequencies of 
two orthogonal cryogenic optical resonators subject to Earth's rotation over ~ 1 year. For a possible 
anisotropy of the speed of light c, we obtain Agc/co = (2.6 ± 1.7) • 10 -15 . Within the Robertson- 
Mansouri-Sexl test theory, this implies an isotropy violation parameter /3 — 8— ^ = ( — 2.2±1.5)-10 , 
about three times lower than the best previous result. Within the general extension of the standard 
model of particle physics, we extract limits on 7 parameters at accuracies down to 10 -15 , improving 
the best previous result by about two orders of magnitude. 
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Special relativity (SR) underlies all accepted theories 
of nature at the fundamental level. Therefore, it has been 
and must be tested with ever increasing precision to pro- 
vide a firm basis for its future application. Such tests 
are also motivated by the efforts to unify gravity with 
the other forces of nature, one of the outstanding open 
challenges in modern science. In fact, many currently 
discussed models of quantum gravity do violate the prin- 
ciples of SR. In string theory, for example, violation of 
Lorentz invariance might be caused by spontaneous sym- 
metry breaking 0; in loop gravity, effective low-energy 
equations (e.g., modified Maxwell equations Q) that vi- 
olate Lorentz invariance have been derived. Since the 
natural energy scale for quantum gravity is the Planck 
scale, E p ~ 10 19 GeV, a direct experimental investigation 
of quantum gravity effects is not feasible. It is, however, 
possible to search for residual effects at an attainable en- 
ergy scale in experiments of outstanding precision. High- 
precision bounds on Lorentz violation might give valuable 
hints for or against particular models of quantum gravity. 

A sensitive probe for Lorentz violation (in electrody- 
namics) is the Michelson-Morley (MM) experiment [3j, 
that even predated the formulation of SR. It tests the 
isotropy of the speed of light, one of SR's foundations. 
In the classic setup, one compares the speed of light c 
in two orthogonal interferometer arms by observing the 
interference fringes. If c depends on the direction of prop- 
agation, the fringes move if the setup is rotated (using, 
e.g., Earth's rotation or a turntable). In 1979, Brillet and 
Hall H introduced the modern technique of measuring a 
laser frequency stabilized to a resonance of an optical ref- 
erence cavity. The frequency of such a resonance is given 
by Vc&v — mc/(2L), where L denotes the cavity length 
and m = 1, 2, . . . the mode number. Thus, a violation of 
the isotropy of c can be detected by rotating the cavity 
and looking for a resulting variation of i^av by comparing 
the frequency of the cavity-stabilized laser against a suit- 
able reference. In our experiment (Fig. !lj), we compare 
the frequencies v x and v y of two similar cavities oriented 



in orthogonal directions, in analogy with the classical in- 
terferometer tests. Compared to the single-cavity setup, 
this arrangement doubles the hypothetical signal ampli- 
tude and provides some common-mode rejection of sys- 
tematic effects. 

Our results are obtained making use of the high dimen- 
sional stability of cryogenic optical resonators (COREs). 
Constructed from crystalline sapphire, COREs show a 
low thermal expansion coefficient (10~ 10 /K at 4.2 K) 
and a remarkable absence of creep (i.e., intrinsic length 
changes due to material relaxation); upper limits for 
the CORE frequency drift are < 2 kHz/6 months and 
< O.lHz/h .5], which makes COREs particularly well 
suited for high-precision measurements such as relativity 
tests 0,0- For the MM experiment, it allows us to rely 
solely on Earth's rotation. This avoids the systematic 
effects associated with active rotation, which previous 
experiments had to use to overcome the creep of room 
temperature resonators made from glass ceramics, e.g. 
ULE (ultra-low-expansion) , on the time-scale of a day. 

A suitable theoretical framework for analyzing tests 
of Lorentz invariance is the general Standard Model 
Extension (SME) d 0. A Lagrangian formulation of 
the standard model is extended by adding all possible 
observer Lorentz scalars that can be formed from known 
particles and Lorentz tensors. In the Maxwell sector, 
the Lagrangian is C = \F av F^ + \{k F ) KXal/ F^F^ 
[lfj. The tensor (kp) K Xfj,u (the greek indices run from 
0...3) has 19 independent components; they vanish, if 
SR is valid. 10 of its components, that can be arranged 
into traceless symmetric 3x3 matrices {h e +) AB and 
(k -) AB , describe polarization-dependent effects. They 
are restricted to < 2 • 10~ 32 by polarization measure- 
ments on light from astronomical sources || and are 
assumed to be zero in the following. The remaining 
9 components describe boost invariance and isotropy 
of c and can be arranged into traceless 3x3 matrices 
(k e ^) AB (symmetric) and (k a+ ) AB (antisymmetric) plus 
one additional parameter. They lead to a shift 8v of the 
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resonance frequency of a cavity 9] with a characteristic 
time signature and can thus be measured in cavity 
experiments. 

For our experiment (Fig. 0), we use two L = 3 cm long 
COREs that feature linewidths of 100 kHz and 50 kHz, 
respectively, (Finesses ~ 10 5 ), located inside a liquid he- 
lium (LHe) cryostat with a liquid nitrogen (LN2) shield. 
Refills and the evaporation of coolants cause mechanical 
deformations of the cryostat, which change the resonator 
position. An automatic beam positioning system actively 
compensates for these movements (Fig. [y. The frequen- 
cies vi as of two diode-pumped Nd:YAG lasers at 1064 nm 
are stabilized ("locked") to resonances of the COREs us- 
ing the Pound-Drever-Hall method. The Nd:YAG laser 
crystal strain (generated by a piezo attached to the crys- 
tal) and temperature were used for tuning the frequency. 
The phase modulation at frequencies f m around 500 kHz 
is also generated by crystal strain modulation using me- 
chanical resonances of the piezo. The light reflected from 
the COREs is detected inside the cryostat (Fig. ^ on 
Epitaxx 2000 InGaAs detectors. Down-conversion to DC 
at 3/ m rather than l/ m reduces the influence of resid- 
ual amplitude modulation and gives a higher signal-to- 
noise (S/N) ratio for the same circulating laser power 
inside the cavity. For the detector signal, amplifiers con- 
sisting of 8 paralleled BF1009 dual-gate MOSFETs pro- 
vide low current noise in spite of the high detector ca- 
pacitance (500 pF). At ~ 80 nW laser power impinging 
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FIG. 1: Setup. Inside a 4K cryostat, two COREs are lo- 
cated in a copper block to provide common-mode rejection 
of thermal effects. LN2 is refilled automatically every ~ 3h, 
LHe manually every ~ 2 days. Laser beams are coupled to 
the COREs via windows, with polarizers P and lock detectors 
PD inside the cryostat. For active beam positioning, beams 
pass through galvanomter (G) mounted glass plates. The hor- 
izontal and vertical displacements are adjusted to maximize 
coupling into the cavities, as measured using the 2/ m signal 
from the detector in reflection. Not shown is the setup for 
measuring the frequency difference v m — v v of the lasers, in 
which the beams are overlapped on a high-speed photode- 
tector and the beat frequency is measured against a quartz 
oscillator stabilized to the GPS. 



on the COREs, we achieve an error signal S/N ratio of 
~ 1.5 • 10 4 • y/rjs. Between 80 . . . 200 nW have been used, 
minimizing the change of z/ cav due to laser heating of the 
COREs (<~ 10Hz///W). The error signal is generated us- 
ing mini-circuits ZHL-32 high-level amplifiers and SAY-1 
23 dBm double-balanced- mixers operating at ~ dBm 
RF amplitude and thus well below saturation. This pro- 
vides highly linear operation and proved very important 
for a low sensitivity to systematic disturbances. On a 
time scale of minutes, we reach a minimum relative fre- 
quency instability of the lasers locked to the COREs of 
7 ■ 10" 16 , referred to a single CORE (Fig. 01. Such a 
level is reached by the best ULE-cavity stabilized lasers 
flll | only if a large linear drift ~ 2Hz/s is subtracted. 

A number of systematic disturbances (e.g. residual 
amplitude modulation, parasitic etalons, or mixer offset 
voltages) that cause lock-point shifts are compensated 
for using a new technique: By an additional phase mod- 
ulation of the laser beams, a part of the laser power is 
shifted into side-bands, thus reducing the amplitude of 
the useful part of the error signal. If the parasitic ef- 
fects do not have sufficient frequency selectivity to dis- 
criminate these sidebands against the carrier, modulat- 
ing the amplitude of the sidebands makes the lock-point 
shifts time-dependent, so they can be identified and com- 
pensated for [l2j . This proved instrumental in achieving 
the reliable and repeatable laser system performance re- 
quired for accumulating enough data. 

Except for a 10 day break around New Year 2002, 
the COREs were operated continuously at 4.2 K over 
more than one year. Usable data (discounting data sets 
shorter than 12 hours and data taken during adjustment 
or LHe refills) starts June 19, 2001 and was acquired over 
390 days until July 13, 2002. A total of 146 data sets 
of 12 h to 109 h in length, totaling 3461 h, are available 
(Fig. |21 and Fig. 49 almost equally distributed data 
sets are longer then 24 h. For extracting results, simul- 
taneous least-squares fits with a constant offset, a linear 
drift, and the amplitude of a sinusoidal signal at fixed 
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FIG. 2: Left: Typical data set fitted with a 12 h sinewave 
amplitude, a linear drift, and a constant offset. Peaks occur 
every few hours due to automatic LN2 refills. Right: Root 
Allan variance calculated from this data (upper curve), and 
from a quiet part between two LN2-refills (118 minutes start- 
ing at 555.87 days; lower curve). 
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frequency and phase as suggested by the test theory are 
performed. Before fitting, the data sets are divided into 
subsets of 12 h (or 24 h for the 24 h signals). This drops 
a fraction of the data, but makes the resulting sinewave 
amplitudes independent of offsets in the data. We obtain 
199 fits of 12 h data subsets (Fig. |3J). 

The individual fit results (Fig. |3J) are combined for the 
final result by coherent (vector) averaging. In the SME, 
the hypothetical Lorentz violation signal (y x — v y )/D — 
J2iAf sinwiT e + Af cosc^T®, where v = {v x + u y )/2, 
has Fourier components at 6 frequencies u>i. As defined 
in 0, T e = on March 20, 2001, 11:31 UT. The sig- 
nal components as given in Tab. [flare derived [l3T | from 
Eq. (39) in (the cavity length L is not substantially 
affected by the hypothetical Lorentz violation Be- 
cause our data extends over more than one year, in the 
vector average we can resolve the 6 signal frequencies |l5j 
and extract elements of the traceless symmetric matrix 

a 1.7 ±2.6 -6.3 ±12.4 
1.7 ±2.6 b 3.6 ±9.0 | • 10~ u ' 
-6.3 ±12.4 3.6 ±9.0 -(a + 6) 

with a — b = 8.9 ± 4.9 (in the sun-centered celestial equa- 
torial reference frame of 0). Likewise, 

14 ±14 -1.2 ±2.6 

/,-,,.., = | -14 ±14 0.1 ±2.7 1 1.0 

1.2 ±2.6 -0.1 ±2.7 

for the antisymmetric matrix. One sigma errors are 
quoted. The limits on k a+ are weaker because its ele- 
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ments enter the experiment suppressed by /3© ~ 10~ 4 , 
Earth's orbital velocity. All elements of k q+ and all but 
one element of K e _ are obtained. Compared to |l6|. we 
improve the accuracy by about two orders of magnitude. 
Moreover, while |16| give limits on linear combinations 
of the elements of k a+ , the present experiment allows in- 
dividual determination, because of the > 1 year span of 
our data. 

For comparison to previous work (e.g., 0,|EE10]), 
we also analyze our experiment within the Robertson- 
Mansouri-Sexl test framework 0. It assumes gener- 
alized Lorentz transformations that contain parameters 
a, (3, and S. In a preferred frame £ (usually the cosmic 
microwave background), the speed of light co=const. In 
a frame S moving with the velocity v with respect to E, 
c/co = 1 — (A + Bsin 2 9)v 2 /cq, where 9 is the angle be- 
tween the direction of c and v. Both A = a — /3 + 1 and 
B — (3 — 5 — i vanish in SR. In this formalism, SR follows 
from (i) Kennedy-Thorndike (KT)-, (ii) MM-, and (iii) 
Doppler shift experiments. The latter determine the time 
dilation coefficient a (= — h in SR) to |a+ || < 8 • 10 -7 
|l8| . KT experiments test velocity invariance of c, usin 
the periodic modulation of v provided by Earth's orbit ^_ t 
or rotation Measuring the frequency of a cryogenic 
microwave cavity (that is proportional to c) against an 
H-maser, obtained A = (-3.1 ± 6.9) • 10" 7 . 

While all three experiments are required for a complete 
verification of SR within this framework, the MM exper- 
iment currently offers the highest resolution. In case of 
a violation of isotropy, B ^ 0, for our experiment we 
obtain a periodic change of the beat frequency v x — v y 
at 2w0 |20j. For such a signal, the experiment yields an 
amplitude 1.03 ±0.53 Hz, or B = (-3.1 ±1.6) • 10~ 9 . As 
the quality of the data is not uniform (Fig. |3J), ta king 
a weighted vector average is more appropriate here |15| . 
Wc divide the data into the intervals A-E (Fig|3Jl with 
approximately uniform data quality within each. The 
averages over the intervals are then combined to the fi- 
nal result, weighted according to their standard error. 
This leads to a signal of 0.73 ± 0.48 Hz, or 



B = 8-6-- 
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-2.2 ± 1.5) • 10 



-9 



which we regard as the final result within the RMS 
framework. It has an inaccuracy about three times lower 
than the best previous result B = (3.0 ±4.9) • 10~ 9 [ll2j]|. 



FIG. 3: Fit results (upper part; dots) and fit residuals y \ 2 
(lower part; crosses) for the amplitude of the cosine compo- 
nent of an assumed isotropy violation signal at 2lj§ . Strongly 
perturbed data sets are omitted when ^p^? > 20 Hz (25 sub- 
sets). In 8 of these good fit was achieved after dis- 
carding the first 1-3 hours (which likely contain residual per- 
turbations from a preceding LHe refill); these were then re- 
introduced to the analysis. The histograms (bin size =1 Hz) 
on the right hand side give the distribution of the fit results 
and residuals. 



In summary, we performed a modern Michelson- 
Morley experiment by comparing the frequencies of two 
crossed cryogenic optical resonators subject to Earth's 
rotation over a period spanning more than one year. 
Within the Robertson-Mansouri-Sexl framework, our 
limit on the isotropy violation parameter is about three 
times lower than that from the classic experiment of Bril- 
let and Hall 4] . Moreover, we obtain limits on seven pa- 
rameters from the comprehensive extension of the stan- 
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TABLE I: Signal amplitudes A?, A? from the SME, based on 0, and fit results. w ffi « 27r/(23h56min) and fig, = 2n/l year 
denote the angular frequencies of Earth's sidereal rotation and orbit; x ~ 42.3° is the colatitude of Konstanz, r\ « 23.4° the 
angle between the ecliptic and Earth's equatorial plane. The signal includes contributions of order 1, or suppressed by either 
the Earth's orbital velocity /3® ~ 10 -4 or by the velocity of the laboratory due to Earth's rotation (3l ~ 10~ 6 . For each U)i, 
we include only the largest term, in effect dropping all terms proportional to /3l- A fit result that has been used to extract a 
parameter of the SME is set in bold face. The unused fit results lead to additional (but weaker) limits on the elements of k +- 
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dard model [alii, down to 10~ 15 , about two orders of 
magnitude lower than the only previous result |l6| . 
While the high long-term stability of COREs allows one 
to use solely Earth's rotation, active rotation could still 
improve the accuracy significantly. Due to the low drift of 
COREs compared to ULE cavities, the optimum rotation 
rate would be relatively slow, which is desirable for mini- 
mizing the systematics. At a rate of ~ 0.2/min (~ 5/min 
was used in [4|), one could utilize the optimum ~ 7- 1CP 16 
frequency stability (Fig. EJ of the COREs, more than 10 
times better than on the 12 h time scale used so far. Since 
~ 500 measurements (2 per turn) could be accumulated 
per day, one should thus be able to reach the 10 -17 level 
of accuracy. Further improvements include fiber coupling 
and COREs of higher finesse. Finally, space based mis- 
sions with resonators are currently studied (OPTIS 
and SUMO 0). 
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